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We report the first beyond-mean-field study of low-lying parity-doublet states in 224 Ra by extend¬ 
ing the multireference relativistic energy density functional method to include dynamical correlations 
related to symmetry restoration and quadrupole-octupole shape fluctuation with a generator coor¬ 
dinate method combined with parity, particle-number, and angular-momentum projections. We 
clarify full microscopically that the origin of spin-dependent parity splitting in low-spin states is 
related to the octupole shape stabilization of positive-parity states, the dominated shapes of which 
drift gradually to that of negative-parity ones. 

PACS numbers: 21.10.-k, 21.60.Jz, 21.10.Re 


The existence of octupole shaped nuclei in intrin¬ 
sic frame has been suggested very early from the ob¬ 
served low-lying parity doublets connected with strong 
electric odd-multipole (EX) transition strengths [1, 2], 
Among them, 224 Ra was recently suggested to be a stable 
pear shaped nucleus based on the measured EX transi¬ 
tions [3]. Nevertheless, as pointed out early [4] and shown 
in this paper that the effects stemming from static oc¬ 
tupole shapes on EX transitions can be almost equally 
well described by including dynamical octupole defor¬ 
mations around the equilibrium shape. Moreover, the 
spectroscopy of reflection-asymmetric diatomic molecules 
suggests the existence of a rotational band with alternat¬ 
ing parity, which has indeed been found in some candi¬ 
dates of octupole shaped nuclei at high-spin states. How¬ 
ever, an evident energy splitting is observed in the odd 
and even parity states of these nuclei at low-spin region. 
In particular, this parity splitting is gradually increased 
with the decrease of spin. This peculiar feature indicates 
the possible existence of large quantum shape fluctua¬ 
tion in octupole shapes at low-spin states. Therefore, 
whether atomic nuclei possess dynamical or stable oc¬ 
tupole shapes in low-spin states remains an open ques¬ 
tion, which is relevant not only for deepening our un¬ 
derstanding on nuclear structure at low energy, but also 
in searching for CP-violating Schiff moment which is ex¬ 
pected to be amplified in octupole deformed nuclei [5]. 

Numerous studies have been performed to search 
for octupole deformed nuclei with selfconsistent mean- 
field approaches [6-12], and to understand nuclear col¬ 
lective octupole excitations with energy density func¬ 
tional (EDF) based collective Hamiltonian [13-16], EDF 
mapped interacting boson model [17], and generator co¬ 
ordinate method (GCM) [18-20]. These studies have 
greatly deepened our understanding on the structure of 
octupole candidate nuclei. However, the onset of unusual 
spin-dependent parity splitting has not been clarified. 


The combination of energy density functional approach 
with GCM provides a powerful tool to study nuclear 
shape fluctuations. Unfortunately, angular-momentum 
projection has not yet been implemented in this frame¬ 
work for octupole nuclei and the stabilization of nuclear 
octupole shape against rotation remains to be examined. 

In the meantime, cranked mean-field approaches have 
also been adopted to study nuclear octupole deforma¬ 
tion in different spin states [21, 22]. A spin-dependent 
staggering in deformations and the energy splitting of 
odd and even parity states were observed in the N = 
88 isotones [22], where a phase transition picture from 
octupole vibration to octupole deformation was sug¬ 
gested. However, the effects of shape fluctuation and 
angular momentum conservation remain to be taken 
into account. On the other hand, to explain the spin- 
dependent parity splitting, Jolos et al. proposed a pic¬ 
ture of rotation-induced second-order phase transition 
from dynamical octupole to static octupole shapes based 
on a one-dimensional collective model with phenomeno¬ 
logical spin-dependent potentials [23], which seems to 
be consistent with semiclassical analysis of available 
data [1, 2, 24, 25]. In contrast, Shneidman et al. pro¬ 
posed a dynamical cluster picture with oscillations in 
mass asymmetry coordinate based on a phenomenolog¬ 
ical dinuclear model [26]. In these phenomenological 
studies, however, the EA transition properties were sel¬ 
dom examined and the microscopic foundation of their 
model assumptions was not clear. In view of all these 
facts, it is important and timely to establish a full micro¬ 
scopic approach to study possible octupole shaped nuclei 
with alternating parity states and to unveil the underly¬ 
ing mechanism responsible for the spin-dependent parity 
splitting at low-spin region. 

To this end, we extend the selfconsistent relativistic 
energy density functional method [27-29] to include dy¬ 
namical correlations related to symmetry restoration and 
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quadrupole-octupole shape fluctuation with the GCM 
combined with parity, particle-number, and angular- 
momentum projections. The symmetry conserved wave 
function for the low-lying collective states is constructed 
by superposing a set of quantum-number projected 
nonorthogonal mean-field reference states | q) around the 
equilibrium shape 
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where the generator coordinate q stands for the dis- 
critized deformation parameters {fa, fa} of the reference 
states from deformation constrained selfconsistent mean- 
field calculation based on a universal relativistic energy 
functional PC-PK1 [30]. The P G s (G = J,ir,N,Z) are 
projection operators [31]. The weight function is 

determined by the Hill-Wheeler-Griffin equation, 

E [*<?«» - fir =° ( 2 ) 

K b 

where the hamiltonian kernel and norm kernel 

zr K J * Kb are given by, 

= (qa\6P J KaKb P N P z P*\q b ) (3) 

with the operator O representing H and 1, respectively. 

In the present work, the mean-field wave functions |g) 
are obtained by selfconsistent deformation constrained 
calculation, where the Dirac equation of single-particle 
wave function is solved in a three dimensional isotropic 
harmonic-oscillator basis with fourteen major shell. Al¬ 
though this basis space is not sufficient for the absolute 
energies of states, it provides reasonable convergent solu¬ 
tions to nuclear spectroscopic properties, including exci¬ 
tation energies and transition strengths. Axial and time- 
reversal symmetries are imposed to reduce the compu¬ 
tational burden. Moreover, since the BCS method and 
Bogoliubov transformation give similar good description 
for the low-lying states of nuclei not far from stability 
line if the pairing strengths are chosen properly [32], we 
adopt the computation less demanded BCS method to 
take into account the pairing correlations between nu¬ 
cleons with the strength parameters of zero-range pair¬ 
ing forces chosen according to the PC-PK1 parametriza- 
tion [30]. The PfafEan techniques [33, 34] are imple¬ 
mented to calculate norm overlaps in the norm kernel. 
The hamiltonian kernel is calculated with mixed-density 
prescription [35]. The numbers of mesh points in rota¬ 
tion Euler angle and gauge angle are chosen as 16 and 
9, respectively, in the angular momentum and particle 
number projections. After convergence check, we finally 
choose 39 reference states in the (fa, fa > 0) deformation 
plane in the GCM calculation. The configurations with 
fa < 0 are included automatically by using the parity 
operator. 


FIG. 1: (Color online) (a) Mean-field energy surface of 224 Ra 
in fa-fa deformation plane (normalized to energy minimum), 
where two neighboring lines are separated by 0.5 MeV. (b) 
Total nucleon density of the energy-minimum state with fa = 
0.179, fa = 0.125, fa = 0.146 in y = 0 plane. 

It is worth mentioning that the GCM based on non- 
relativistic Skyrme [36] or Gogny [37] forces has recently 
been extended for odd-particle states [38] and cranked 
states [39], respectively with great efforts. However, 
these approaches are presently limited to quadrupole 
shapes. The present method is a further extension of 
our previous beyond mean-field method for triaxial nu¬ 
clei [40] by allowing for octupole shape degree of freedom 
and provides the most advanced beyond mean-field cal¬ 
culation for nuclear octupole excitations. 

Figure 1(a) displays mean-field energy surface of 224 Ra 
in fa-fa deformation plane, where the deformation pa¬ 
rameters fa are defined as 

fa = 3^x<5k A ^ok), R = 1.2 A 1 / 3 , (4) 

with A being mass number of the nucleus. The topog¬ 
raphy of the mean-field energy surface in Fig. 1 is essen¬ 
tially the same as that obtained in Ref. [17] with a rel¬ 
ativistic Hartree-Bogoliubov (RHB) model. The global 
energy minimum is located at the configuration with de¬ 
formations fa = 0.179, fa = 0.125 ,/?4 = 0.146, which 
are slightly larger than the values fa = 0.154, fa = 
0.097, /?4 = 0.15 in Ref. [3]. Total nucleon distribu¬ 
tion of this energy-minimum configuration is depicted 
in Fig. 1(b), exhibiting a well-developed pear shape. 
The energy difference between the energy minimum and 
the lowest reflection-symmetric configuration is, how¬ 
ever, not large enough to prevent shape fluctuation. The 
present GCM calculation shows a broad distribution of 
the collective wave function in fa-fa plane for ground 
state, consistent with the finding in Ref. [20]. With 
the increase of spin, the spread collective wave function 
of positive-parity states becomes gradually concentrated 
around the energy minimum and close to that of negative- 
parity states, which is more stable against rotation. 

The calculated low-lying energy spectra are compared 
with available data in Fig. 2. The spectra and the E2 and 
E3 transitions can be reproduced reasonably well using 
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FIG. 2: (Color online) Low-lying energy spectra for 224 Ra. 
The available data from Ref. [3] and the results from full- 
configuration and the single energy-minimum configuration 
calculations are shown in (a), (b) and (c), respectively. The 
numbers on arrows are E2 (red color) and E3 (blue color) 
transition strengths (Weisskopf units). 


only the energy-minimum configuration (cf. Fig. 2(c)). 
However, the predicted alternating parity rotation band 
as expected for a stable octupole shaped nucleus is not 
supported by the data. The energy displacement be¬ 
tween the parity doublets can only be reproduced in the 
GCM calculation (cf. Fig. 2(b)) by mixing the con¬ 
figurations around the equilibrium shape, which does 
not affect appreciably the good description for the EA 
strengths. As clearly shown in Fig. 3, the gradually re¬ 
duced odd-even staggering Rj/ 2 with increasing J is asso¬ 
ciated with the similar behavior of average deformation 
parameters (f3\) = Y^ q I /5a 1 1 77 1 2 ; where the orthonor¬ 
mal collective wave function is constructed as [31] 
9 J c?((la) = J2q b [^ j7I ] 1 q / a 2 qb fq™- The octupole defor¬ 
mations of negative-parity states are almost constant 
(around 0.144) with increasing spin J, cf. Fig. 3(c), 
while those of positive-parity states increase smoothly 
from </? 3 ) = 0.101 (J = 0) to (/? 3 ) = 0.141 (J = 10 h). A 
similar but much less evident staggering is shown in (/3 2 ) 
varying in between 0.180 and 0.191. 

It is worth mentioning that in contrast to the predic¬ 
tions based on the energy-minimum configuration, the 
full configuration-mixing calculation predicts a steady in¬ 
crease (but also with somewhat odd-even staggering) of 
electric octupole transition moments Q 3 (J —> J — 1) and 
Q 3 (J —1 J — 3). However, the available data on the 
EA transitions in 224 Ra is limited and can hardly distin¬ 
guish these two situations. Moreover, the ratio of electric 
dipole moment Qi(J —> J — 1) to quadrupole moment 
Q 2 (J —> J — 2) is similar in both calculations and not 
much dependent on spin J, which is consistent with the 
available data [24], In short, Fig. 3 presents us a novel 
spin-controlled shape stabilization picture that rotation 
brings the dominate shape of positive-parity (even-,/) 



FIG. 3: (Color online) (a) Excitation energy ratio Rj /2 = 
E X (J)/E x ( 2 + ), (b) normalized staggering amplitude Sj / 2 = 

\ E *(J)-4rhE x (J-l)-^ T iE x {J + l)\/E x (2+) [25] and (c) 
average deformation (/ 3\) with A = 2, 3 as functions of spin J 
in 224 Ra. 


states close to that of negative-parity (odd-J) states. A 
further high-precise measurement on the Q 3 (J —> J — 1) 
transition moments should be helpful to confirm this pic¬ 
ture derived from the excitation energies. 

To understand the connection between spin-dependent 
parity splitting and shape evolution in a simple and intu¬ 
itive way, we carry out the GCM calculations by mixing 
the octupole configurations at either /3 2 = 0.10 or 0.15. 
Fig. 4(a) displays the corresponding excitation energy ra¬ 
tio Rj /2 and normalized staggering amplitude Sj/ 2 ■ The 
gradually reduced energy staggering is shown in both 
cases. The mechanism responsible for this behavior is 
exhibited in Fig. 4(c) and (d) for the projected energy 
surfaces and in Fig. 4(e) and (f) for the collective wave 
functions of each spin-parity states. We note that the 
peak at (3 3 = 0.02 in Fig. 4(d) is due to the collapse of 
pairing correlation between protons and this configura¬ 
tion turns out to have a minor influence on the GCM 
results. 

In the projected GCM calculation, energy splitting be¬ 
tween AJ = 1 states comes mainly from the dynamical 
correlation energies (DCE) from angular-momentum and 
parity projections. The DCE from parity projection is 
decreasing to zero with increasing / 3 3 as the norm overlap 
(p 2 p 3 \P\p 2 p 3 ) —> 0, cf. Fig. 4(b). Moreover, the larger 
the /3 2 is, the faster the norm overlap approaches zero 
with f3 3 . Therefore, the energy splitting between A J = 1 
states decreases with increasing of both j3 2 and / 3 3 , as 
shown in the projected energy surfaces, cf. Fig. 4(c) 
and (d). An alternating-parity projected energy surfaces 
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FIG. 4: (Color online) (a) Excitation energy ratio Rj /2 
and (inset) staggering amplitude Sj/ 2 ; (b) norm overlap 
(/? 2 /? 3 1|/? 2 /? 3 ) of parity operator in mean-field states; (c) and 
(d) mean-field and projected energy curves onto particle num¬ 
bers N, Z , angular momentum J, and parity n = (—1) J (solid 
lines for even J, dashed lines for odd J). The mean-field (or 
projected) energy curves are normalized to the 0 + state at 
fa = 0 (or minimum). The mean-field energy surfaces in (c) 
and (d) are shifted down by 3.43 MeV and 4.15 MeV, respec¬ 
tively. (e) and (f) absolute value of collective wave functions 
for parity-doublet states. Those for odd-7r ( 1 _ , 3 “ ,..., 9 - ) 
states are close to each other. 


is shown in sufficient large deformation regions, where 
the energy splitting between A J = 1 states is originated 
almost purely from rotation, as expected for stable oc- 
tupole shaped nuclei. In contrast, for the configurations 
in relatively small deformation regions, such as fa < 0.05, 
the nonzero overlap (fafa\P\fafa) leads to a large energy 
splitting between A J = 1 parity doublets and thus to a 
large staggering amplitude. 

In particular, Fig. 4(c) and (d) shows that the pro¬ 
jected energy surfaces with positive parity evolve evi¬ 
dently with increasing spin, presenting a transition from 
weakly deformed octupole shape to a large deformed 
one, which is consistent with the behavior of the spin- 
dependent collective potentials introduced phenomeno¬ 
logically in Ref. [23]. This spin-controlled shape stabi¬ 
lization picture is demonstrated more clearly in the dis¬ 
tributions of collective wave functions in Fig. 4(e) and 
(f). With the increase of spin, the dominated configu¬ 
ration of positive-parity state drifts gradually from the 
weak octupole configurations {fa ~ 0.05) to those with 


large octupole shapes {fa £ [0.10,0.15]). This behav¬ 
ior is probably related to the deformation-dependent mo¬ 
ment of inertia, which turns out to increase with fa and 
fa, presenting a peak region around the minimum of the 
mean-field energy surface in Fig. 1. In contrast, due to 
antisymmetry requirement {fa —» — fa) for the negative- 
parity states, their collective wave functions are zero at 
fa = 0, concentrated to a large octupole deformed con¬ 
figurations {fa £ [0.10,0.15]) and are hardly affected by 
rotation. It is coincident with the evolution picture of the 
collective wave functions with spin from the full GCM 
calculation. 


In conclusion, we have reported the first symmetry- 
conserved beyond-mean-field study of low-lying parity 
doublets states in 224 Ra with a state-of-the-art mul¬ 
tireference relativistic energy density functional method, 
where the dynamical correlations related to restoration 
of broken symmetries and to fluctuations of quadrupole- 
octupole shapes have been taken into account with 
the exact generator coordinate method combined with 
particle-number, angular-momentum, and parity projec¬ 
tions. Both the energy spectrum and the E2, E3 tran¬ 
sitions are reproduced in the configuration-mixing cal¬ 
culation without introducing any phenomenological pa¬ 
rameter. We have demonstrated full microscopically 
that 224 Ra has dynamical octupole shapes in the low- 
spin positive-parity states and acquires a stable octupole 
shape close to that of negative-parity states after rota¬ 
tion with spin up to ~ 8 h. This novel picture of rotation- 
induced octupole shape stabilization provides a natural 
explanation for the spin-dependent parity splitting in the 
excitation energies and is expected to be a common phe¬ 
nomenon in some actinides and rare-earth nuclei. Fi¬ 
nally, we point out that with our microscopic results as 
the inputs of the coupled-channel calculation for nuclear 
fusion [41], one may observe some interesting differences 
on the barrier distribution compared with that based on 
rotational or vibrational limits [42]. Moreover, octupole 
deformed odd-mass nuclei have enhanced time-reversal 
violating nuclear Schiff moments and thus are relevant 
for measuring atomic electric-dipole moments [43]. The 
present study provides a starting point to examine the 
dynamical correlation effects on nuclear Schiff moments. 
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